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Three-dimensional Two-Layer Outer Gap Model: Fermi 
Energy Dependent Light Curves of the Vela Pulsar 
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ABSTRACT 

We extend the two-dimensional two-layer outer gap model to a three-dimensional 
geometry and use it to study the high-energy emission of the Vela pulsar. In this 
model, the outer gap is divided into two parts, i.e. the main acceleration region on 
the top of last-open field lines and the screening region around the upper boundary 
of the gap. In the main acceleration region, the charge density is much lower than 
the Goldreich- Julian charge density and the charged particles are accelerated by the 
electric field along the magnetic field to emit multi-GeV photons. In the screening 
region, the charge density is larger than the Goldreich- Julian value to close the gap 
and particles in this region are responsible for multi-lOOMcV photon emission. We 
apply this three dimensional two-layer model to the Vela pulsar and compare the 
model light curves, the phase-averaged spectrum and the phase-resolved spectra with 
the recent Fermi observations, which also reveals the existence of the third peak 
between two main peaks. The phase position of the third peak moves with the photon 
energy, which cannot be explained by the geometry of magnetic field structure and 
the caustic effects of the photon propagation. We suggest that the existence of the 
third peak and its energy dependent movement results from the azimuthal structure 
of the outer gap. 
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1 INTRODUCTION 

The Vela pulsar is one of the few 7-ray pulsars, whose 7-ray 
radiations are so intense that could be detected before the 
launch of Fermi 7-ray telescopes. For this reason, the ob- 
served spectra and the pulse profiles of the Vela pulsar have 
been studied extensively before the launch of Fermi (Ro- 
mani, 1996; Dyks, Harding & Rudak, 2004; Takata, Chang 
& Shibata 2008). Its 7-ray spectrum between lOOMeV to 
several GeV is difficult to be explained in terms of simple 
curvature radiation. Takata, Chang & Shibata (2008) pro- 
posed that the component around lOOMeV of the Vela pul- 
sar, as well as Geminga (Takata & Chang, 2009), is due to 
the synchrotron radiation of the incoming particles in the 
outer gap. By solving the evolution of the Lorentz factor 
and the pitch angle in the two-dimensional electrodynamics 
model (Takata, Shibata & Hirotani, 2004), they found that 
the synchrotron radiation of the incoming particles dom- 
inates the radiation around 10 MeV to lOOMeV and the 
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curvature radiation of the outgoing particles dominates the 
radiation around lGeV. 

Zhang & Cheng (1997) proposed a self-sustained thick 
outer gap model of 7-ray emission from the rotation-powered 
pulsar. They pointed out that the primary e ± pairs in the 
steady state have a power-law distribution and used this 
distribution to calculate the spectra of some 7-ray pulsars, 
including Geminga-like and Vela-like pulsars. Also based on 
this distribution of primary pairs, Zhang & Cheng (2001) 
applied a three-dimensional pulsar magnetosphere model to 
explain the high-energy emission from the Geminga pulsar 
with a thick outer gap. In this calculation, the high-energy 7- 
rays are produced by the accelerated particles with a power- 
law energy distribution via curvature radiation inside the 
outer gap. 

After the launch of Fermi LAT, the bottleneck of the 
research of high energy radiation from pulsar is broken. Just 
in one year, this high-quality 7-ray telescope has measured 
46 7-ray pulsars (Abdo et al. 2010a). We developed a two- 
layer model (Wang et al., 2010) to study the phase averaged 
spectra of the mature pulsars, including the Vela pulsar, in 
the first catalogue of 7-ray pulsars of Fermi. In this simple 
two-dimensional model, the outer gap extending outwards 
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from the null charge surface is considered as a superposi- 
tion of a main acceleration region and a screening region. 
In the main acceleration region, which is close to vacuum, 
the electric field is strong, and the photons generated in this 
region dominate the radiation of around GeV range. In the 
screening region, the large numbers of the pairs produced by 
the pair-creation process stop the growth of the acceleration 
region and they produce the photons with around lOOMeV 
range. Although this two-dimensional two-layer model could 
explain the phase-average 7-ray spectra of the mature pul- 
sars, it cannot be used to study the phase-resolved spec- 
tra and the light curves. Furthermore the high quality data 
of Fermi has clearly shown the existence of a third peak 
(bump) at the trailing part of the first peak, which was ex- 
plained as an effect of geometry of the magnetic field (Dyks 
& Rudak, 2003; Fang & Zhang, 2010). However, the ob- 
served result provided by Fermi (Abdo et al, 2010b) shows 
that, as the energy increases, this third peak shifts towards 
the second peak. This indicates that unlike those two main 
peaks this third peak does not result from the caustic effect 
of photon propagation and/or the structure of the magnetic 
field. 

The high-energy radiation processes of pulsars have 
been studied based on polar cap models (Ruderman & 
Sutherland 1975; Daugherty & Harding, 1982), the slot 
gap models (Arons 1983; Harding, Usov & Muslimov 2005; 
Harding et al. 2008) and the outer gap models (Cheng, 
Ho & Ruderman 1986a,b; Romani 1996; Hirotani 2008; 
Takata & Chang 2009). The different acceleration models 
have predicted different properties of the 7-ray emission 
properties from the pulsar magnetosphere. To discriminate 
among these existed emission models, it is better to com- 
pare the predicted phase resolved spectra with the observed 
ones, which contains the most detailed information about 
the structure of pulsar magnetosphere and the acceleration 
mechanism. 

In this paper, we extend our study of two-layer model 
from two-dimensional to three-dimensional, which includes 
the effects of the inclination angle, viewing angle and the az- 
imuthal structure of the gap. This three-dimensional model 
can provide the phase-resolved spectra and the energy- 
dependent light curves that cannot be given by the two- 
dimensional one. We will use it to study the observed en- 
ergy dependent light curves of the Vela pulsar measured by 
Fermi LAT (Abdo et al, 2010b). 

In section 2, we describe our theoretical model, which 
contains the electrodynamics of the outer gap, the emission 
geometry and the method to calculate the spectra and the 
light curves. In section 3, we will use our three-dimensional 
model to explain the energy dependent light curves of the 
Vela pulsar and discuss the conditions of the existence of and 
the movement of the third peak as well as other observed 
spectral properties. We gives a brief summary in section 4. 



2 THEORETICAL MODEL 

2.1 Three dimensional gap structure 

We discuss the 7-ray emissions from the outer gap accelera- 
tor in the pulsar magnetosphere, and assume that a strong 
acceleration region is extending above the last-open filed 
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Figure 1. The simplified distribution of the charge density (solid 
line) and the corresponding accelerating potential (dot-dashed 
line) of the two-layer outer gap. 



lines and between the null charge surface of the Goldreich- 
Julian charge density and the light cylinder. Although the 
dynamical studies (e.g. Takata, Shibata & Hirotani, 2004) 
predict that the inner boundary of the active gap can extend 
inwards from the null charge surface, most of the power of 
the gap are released beyond the null charge surface. Because 
the pairs are always created on convex side of the magnetic 
field lines, the charge number density increases as the height 
measured from the last-open field lines. In the trans-field 
direction of the magnetic field, therefore, as shown in Fig- 
ure [T] the outer gap can be divided into two parts: 1) the 
main acceleration region at the lower part of the outer gap, 
where the charge density is assumed to be ~ 10 % of the 
Goldreich- Julian value and a strong electric field is accel- 
erating the particles to emit GeV photons via the curva- 
ture process, and 2) the screening region around the upper 
boundary, where the growth of the main acceleration region 
in the trans-field direction is stopped by the pair-creation 
process. Such two-layer structure, in the direction along the 
field line, starts from a place between the inner boundary 
and the null charge surface, and ends at the light cylinder. 
Inside the inner boundary and outside the light cylinder, 
there is no accelerating electric field; from the null charge 
surface to the inner boundary, the accelerating electric field 
decreases rapidly to zero. 

Because the number density of the pairs made by the 
pair-creation process increases exponentially in the trans- 
field direction (Cheng, Ho & Ruderman 1986a, b), we use a 
simple step function to approximate the distribution of the 
charge density in the trans-field direction in the poloidal 
plane (the plane where the field lines have same polar angle 
4> P ) as 



p{x,z,<j> v ) = 



pi(x,<pp), if < z s£ h\ (x, <pp) 
p2(x,4> p ), if hi(x,cj> p ) < z < h 2 (x,c/>p) 



(1) 



where x, z and 4>v represent coordinates along the magnetic 
field line, the height measured from the last-open field line, 
and the azimuthal direction. In addition, hi and /12 repre- 
sent the thickness of the main-accretion region and the total 
gap thickness, respectively. As illustrated in Figure [T] in the 
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z-direction, the main acceleration region spans from z = 
to z — hi, where the charge density is substantially smaller 
than the Goldreich- Julian charge density (Goldreich & Ju- 
lian 1969), and the screening region is between hi ^ z ^ h 2 , 
the charge density in it is much higher than pcj in order 
to close the gap. For each poloidal plane (0 p =constant), 
hi(x,<j) p ) as well as h 2 (x,(j>p) is defined with a given mag- 
netic field line in the constant <^-plane so that the ratio 
hi(x,4>p)/h 2 (x,cj> p ) is nearly constant along the field lines. 

Even though the structure of the gap is simplified as de- 
scribed above, the three-dimensional Poisson equation with 
the magnetic field structure is still difficult to be solved 
in comparing with the case of the two-dimensional model 
discussed in Wang et al.(2010). Instead of solving the real 
three-dimensional Poisson equation, we adopt a simplified 
Poisson equation. We calculate the electric field structure at 
a fixed poloidal plane, that is, at each azimuthal coordinate 
by solving the two-dimensional Poisson equation. This ap- 
proximation can be justified if the derivation of the potential 
field in the azimuthal direction is much smaller than that in 
poloidal plane. This procedure may be applicable for the 
thin outer gap accelerator, whose thickness in the poloidal 
plane is much smaller than the width in the azimuthal direc- 
tion. The detail procedure to obtain the three-dimensional 
structure is described as follows. 

In a fixed poloidal plane, the bottom and top bound- 
aries of the outer gap are the last-open field lines and a 
surface layer with separation h 2 made by the magnetic field 
lines respectively. We divide the polar cap into Nb equal 
divisions in the azimuthal direction. If the gap size of the 
azimuthal direction is much larger than the gap thickness in 
the poloidal plane, the gap structure of each slice with a fixed 
azimuthal angle 4> p may be approximately represented by 
the two-dimensional situation shown in Wang et al. (2010). 

The accelerating electric field in each slice of the gap 
is solved by the method described in the previous two- 
dimensional model(Wang et al., 2010). The potential in the 
gap is described by the Poisson equation, X/ 2 & = —inp. This 
potential can be written as <E> = $' + $0, where $0 is the 
co-rotating potential, which satisfies V 2 ^o = —4irpGJ, and 
the potential $' represents the deviation from co-rotating 
potential. This deviation generates the accelerating electric 
field. By assuming the derivation of the potential field in 
the azimuthal direction is small enough to be neglected com- 
pared with that in the poloidal plane, we express the Poisson 
equation of $' as 



dx 2 dz 2 



$' = — 4n(p - 



PGJ) 



(2) 



Because our two-dimensional study predicts that the outer 
gap of the Vela pulsar is thin, the variation of the Goldreich- 
Julian charge density in the trans-field direction (that is z- 
direction) may be able to be neglected so that the pcj can 
be assumed to be constant in z-direction. 

Moreover, we assume that the derivative of the poten- 
tial field in the trans-field (z) direction is larger than that 
along (x-direction) the magnetic field line. This assumption 
may be applicable for the thin outer gap accelerator of the 
Vela pulsar. By ignoring the derivative of the potential field 
along the magnetic field line, we reduce the two-dimensional 
Poisson equation ([2]) to one-dimensional form as 



d 

— $'(x, z, 4> p ) = -4tt[ P (x, z, 4> p ) - paj(x, 4> p )]. 



(■■>>) 



At the lower boundary, where the accelerating electric 
field just begins, the boundary condition is 



&(x,z = O,0 P ) = 0. 



(4) 



At the upper boundary, which is located on the surface of 
fixed magnetic field line, we impose 



<&'(x,z = h2,4> p ) = and E±(x, z = h 2 , <f> p ) = 0. 



(5) 



Using the above boundary conditions and imposing that $' 
and d& /dz are continuous at the height hi, we obtain the 
solution for the Poisson equation (j3]) as 



&(x,z,cf> p ) 



-2n < 



{pi{x,4> p ) - p G j(x,(j) p )}z 
+Ciz, 

for < z < hi 

{p 2 (x,(j} p ) - p G j(x,4> p )} 

(z 2 -hUx,<f> P )) 
+Di(z — h 2 (x, <j) P )), 
for hi z ^ h 2 



(6) 



where 



Ci = 



Di = 



(pi - pGj)hi{hi - 2h 2 ) - (p 2 - PGj)(hi - h 2 y 



(pG.7 - p2)h\ + (pi - p 2 )h\ 



Because the two boundary conditions are imposed on the 
upper boundary (z — h 2 ), the position of the upper bound- 
ary cannot be chosen arbitrarily. On the other hand, the 
condition E±(x,z = h 2 ,<f) p ) = gives the relation between 
the charge densities (pi, p 2 ) and the gap thickness (hi, h 2 ) 



(p 2 - PGj)h 2 + (pi - p 2 )h\ 



0. 



The present model expects that the charge density 
should proportional to the Goldreich- Julian charge density 
in the screening region (see Wang et al. 2010), and there- 
fore we express p(x, z, 4> p ) — pc,j(x, </> p ) ~ g(z,<f>p)pGj(x, 4> p ), 
where 



g(z,<t>p) 



-9i( 



if ^ z hi 
if hi < z ^ hi 



(8) 



We assume that gi > and g 2 > so that \p\ < \pcj\ in 
the main region and \p\ > \pgj\ for the screening region, 
respectively. Substituting equation ((Sj into equation ([7]), we 
obtain 



1 + 



£1 

92 



(0) 



We approximate that the Goldreich-Julian charge den- 
sity as pgj(x) ~ ~^§^ (Cheng, Ho Ruderman 1986a), 
where Q is the angular frequency of the pulsar. Using 
the relations that d(Bh%)/dx ~ 0, d(z/h 2 ) / dx ~ 0, 
d(hi/h 2 )/dx ~ 0, and ds/dx ~ 0, the accelerating electric 
field, -En = —d<&'/dx is written down as 



E\\(x,z,4>p) 



where 



-giz 2 + C[z, 
Q,B I for s: z s; hi 

g 2 (z 2 -hl) + D'i(z-h 2 ), 

for hi < z ^ h 2 



(10) 
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ffifei(fei - 2h 2 ) + ff2(/ti ~ h 2 ) 2 
h 2 

92hl + (gi + g 2 )h\ 



In the three dimensional rotating magnetic field struc- 
ture, the thickness of the gap h 2 is a function of the distance 
to the star along the field line x. We define the gap fraction 
/ measured on the stellar surface as (Wang et al., 2010), 

h 2 (R s ) 



f 



(11) 



where R a is the stellar radius, and r p (4> p ) is the polar cap 
radius and depends on polar angle in three-dimensional mag- 
netic field geometry. Note that because the electric field E^ 
is proportional to Bh 2 , which is almost constant along the 
field line for the dipole field, it can be found that E\\ oc f 2 . 
This indicates that the strength of the electric field and re- 
sultant emissivity of the curvature radiation increase with 
the fractional gap thickness. 

The height(z) measured from the last-open field lines 
is a function of distance (a;) along the magnetic field line. 
In order to make it easy to distinguish the two layers in 
any x, we introduce the factor a to represent the mag- 
netic field lines at a given layer and take a — 1 for the 
last-open field lines. We first determine the coordinate val- 
ues [Xo((f> p ),Yo((j>p), Zo((f>p)] of the last-open field lines at 
the stellar surface. Then an arbitrary layer of the mag- 
netic field lines with a given "a", whose coordinate val- 
ues [X' {(j>p)y Yo(<j> p ), Z' (<f> p )], can be determined by using 
X' {<t> v ) = aXo(<f> P ), Yi{<t> v ) = aYo^p) and Z' {4> v ) = (R 2 - 
c^X'q — a 2 !^' 2 ) 1 / 2 . The relation between the z and the a is 



approximated as 
i 

Z{x,(jip 



1 



-h 2 (x,4>p) 



(12) 



where a m i n corresponds to the upper boundary of the gap. 
In this paper, the dimensionless thickness a m i„ is assumed 
to be the same for different polar angle (f> p because each layer 
should have a similar polar cap shape on the stellar surface, 
and is chosen as the fitting parameter of the light curve. The 
thickness a m i n determines the width of the pulse peaks and 
the phase separation between two main peaks. 



2.2 Curvature radiation spectrum 

Our calculation method for the phase-average spectrum and 
phase-resolved spectra is based on Tang et al. (2008). In a 
volume element AV, there are AN = AVn particles acceler- 
ated by the electric field that is described by equation (110[) . 
The accelerated particles release the power gained from the 
accelerating electric field, eE^c, through the curvature ra- 
diation process. The total radiation power for each particle 
is l C ur{r) = 2e 2 c7e (r)/3s 2 (r), where f = (x, z, <f> p ). The 
typical local Lorentz factor 7 e of the primary particles can 
be obtained by requiring eE\\c = l cur , which gives 



7 ^ = [ |£!M W]1 /4 



(13) 



The photon spectrum at each position of radiation is 
d 2 N J (E J ,r) _ AN{r)V3e 2 f, 



dE-ydt 



2-kUsE-, 



-F(X), 



(14) 



where % = E 1 /E cur {f), E cur (r) = (3/2)hcrfl(r)/s(r) is the 
characteristic energy of the radiated curvature photons, s(f) 
is curvature radius, and F(x) = J°° K 5 / 3 (^)d^, where K 5 / 3 
is the modified Bessel functions of order 5/3. The local solid 
angle of curvature photon AQ ~ ty9 2 , where 



0(r) = Max{Al/s{r),-i{r)mc 2 /eB{r)s{r)), 



(15) 



where the first term represents the solid angle caused by the 
curvature of the magnetic field line and the second term is 
the minimum emission angle measured from the direction of 
the magnetic field line. Here Al is the gird size of the cell 
along field line. 

With equation JS), we express the number density at 
each point as n(r) — — g(z, 4>p)\- The volume element 

at each position can be calculated 

A,/« _ 2TTr p h(R s )B(R s )Al(x) ^ -v 

^ V V I — N A N B B(r) ' 

where 2irr p h 2 (R 3 )B s is the total magnetic flux through in 
the gap, and Na and Nb are the number of the grids in the 
direction of the trans-field direction (z-direction) and the 
azimuthal direction, respectively. 

The total photon flux received at Earth is 



Ftot(Ej 



1 



E 



d 2 N 1 {E 1 ,f,) 



(17) 



where fi represents the position of the i' cell, from which 
the emission can be observed, and D is the distance to the 
pulsar from Earth. We note that because Ti depends on the 
viewing angle and the inclination angle, the spectrum (|17[) 
is a function of the viewing angle and the inclination angle. 

To calculate the phase-resolved spectra and the light 
curves, the arrival times of the photons are binned by pulse 
phase. For example, Figure [2] represents photon-mapping on 
the plane spanned by the viewing angle and the pulse phase 
using the rotating dipole fields (section l2.3p . For each view- 
ing angle, the number of the photons measured at pulse 
phases between ip± and ip 2 is calculated from 



E-2 



FtotiE^x^ii^^dE. (18) 



2.3 Emission Geometry 

It has been considered that the positions of the first and the 
second peaks of the pulsar are determined by the geometry 
of the magnetic field. In this study, we adopt the rotating 
vacuum dipole field to calculate the light curves and spectra. 
For the rotating dipole, the magnetic field B(r) is given by 

S = f[r-(^ + ^ + A )] _(4 + 4 + A ) (19) 

(Cheng, Ruderman & Zhang, 2000), where ft = 
[i(x sin a cos fit + y sin a sin fit + icosa) is the magnetic 
moment vector, f is the radial unit vector, and ct is the 
inclination angle. The Runge-Kutta method is employed 
to trace out field lines and to find the polar cap rim 
[X {4>p),Y ^p),Z Q {4> P )]. 

Once the upper boundary a m i n is chosen, we trace the 
open field lines defining the outer gap, and calculate the 
emission direction in the observer's frame at each point on 
the field lines. In Tang et al. (2008), the emission direction 
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coincides with the direction of the rotating vacuum dipole 
field defined in the co-rotating frame. In the present three- 
dimensional model, on the other hand, the emission direction 
calculated in only observer's frame with the computation 
method employed by Takata et al. (2007). The curvature 
photons are assumed to be emitted in the direction of the 
particle motion, which can be described as 

v = v p B/B + rxti (20) 

where the first term represents the motion along the mag- 
netic field line, v p is calculated from the condition that 
\v\ = c, and the second term is the drift motion. The polar 
angle to the rotation axis £ of the emission direction and the 
pulse phase ip are calculated from 

f cos £ = ^ , 21 > 

[ tp = - cos^ 1 (v x / v xy ) - f ■ v/R L 

(Yadigaroglu 1997). If a point in a field line satisfies |£ — /3| < 
f(r), the radiation of that point can be seen by the observer, 
where j3 is the viewing angle and <p(r) = Al(r)/s(r). Be- 
cause the positions of the first and second peaks in the light 
curve is also determined by viewing geometry, the predicted 
inclination angle a and the viewing angle /3 are uniquely de- 
termined by comparing the calculated peak separation with 
the observations. 



3 RESULTS AND DISCUSSION 

In this section, we apply our model to the Vela pulsar. In 
our two-layer model, the general properties of 7-ray emis- 
sions from the outer gap are characterized by /, gi and 
hi/ti2- Specifically, the gap fraction / mainly determines 
the cut-off energy in the spectrum. Zhang & Cheng (1997) 
suggested a self-consistent mechanism provided by the X- 
rays emitted by the polar cap that is heated by the return 
current to restrict the gap. Their fractional size of the gap, 
fzc = 0.32P 26 / 21 > predicts that of the Vela pulsar 

is 0.15. In the two-dimensional model (Wang et al. 2010), 
/ = 0.16 was used for fitting the phase-averaged spectrum 
of the Vela pulsar. The current in the main acceleration re- 
gion, 1 — gi, and the ratio between the thicknesses of the 
main acceleration region and that of the whole gap, /11//12, 
determine the photon index of the spectrum together. For 
the Vela pulsar, the fitting of the two-dimensional two-layer 
model gives I - g x = 0.08 and hi /ha = 0.927 (Wang et al. 
2010). 

The position of the upper boundary a m m is treated 
as a model parameter, because it is affected by the actual 
magnetic field structure, which is not understood well. To- 
gether with the inclination angle a and the viewing an- 
gle /?, the position of the upper boundary a m ; n can be 
constrained by the detailed structure of pulsed profile. As 
we show in the following section, we find that the set of 
) = (57°, 80°, 0.935) reproduces well both the 
light curves and the spectra. Figure [2] shows the sky-map of 
the emitted photons from the magnetic surface of a = 0.95 
with a = 57° and /3 = 80° . Although there are many sets of 
the inclination angle and the viewing angle that can provide 
the observed peak separation, our choice of these two angles 
is not arbitrary. The reason for this will be discussed later. 
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Figure 2. The skymap of the radiations outside of the null charge 
surface with inclination angle a = 57°. The x-axis is the pulse 
phase and the j/-axis is the direction of the radiation. The dashed 
line is the viewing angle, which is chosen as /3 = 80° . 

3.1 A Simple Case 

Firstly, we assume constant 1 — gi, hi/h% and / in the az- 
imuthal direction. By fitting the phase-averaged spectrum 
for the Vela pulsar, we obtain 1 — g\ = 0.05, h\/hi = 0.927 
and / = 0.2 for the present three-dimensional model. Here 
we assume the distance of D = 325 pc. 

Figure [3] compares the calculated (histograms) and the 
observes (solid lines) pulse profiles in six different energy 
bands. We can see in Figure [Sjthat the model provided pulse 
profiles are generally consistent with the observations. For 
example, the pulsed profiles in wide energy bands have two 
main peaks with the separation of about 0.42. The positions 
of the two peaks in the light curves are determined by the 
geometry of the magnetic field. 

In the present simple model, however, the calculated 
pulse profile can not explain the existence of the observed 
third peak between two main peaks, as we can see in Fig- 
ure [3] The peaks emerging in the calculated pulse profiles 
of Figure [3] are caused by so called caustic effects (Yadiga- 
loglu and Romani 1995; Cheng, Ruderman and Zhang 2000; 
Dyks, Harding and Rudak 2004), in which more photons 
are observed at narrow width of the rotation phase due to 
the special-relativistic effect (that is, the aberration of the 
emission direction and photon's travel time). In the caus- 
tic model, the phases of the peaks are determined by the 
magnetic field structure and they do not depend on the en- 
ergy band. On the other hand, the observed phase of the 
third peak shifts with the energy bands. The simple caustic 
model can not explain the phase shift of the observed third 
peak. On these ground, we speculate that the existence of 
the third peak and its phase shift is related with more com- 
plex structure of the emission region. 



3.2 More Realistic Case 

In this section, we discuss the third peak with a more com- 
plex gap structure by taking into account the azimuthal 
structure of the fractional gap thickness /, the ratio of the 



6 Y. Wang, J Takata & K.S. Cheng 




0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 



3000 - 






1 0.1-0.3 GeV 


120 - 




0.02-0.1 GeV 


3000 - 


2500 - 








100 - 






2500 - 


2000 - 






80 - 






2000 - 


1500 - 


1 






60 - 








1 500 - 


1000 - 








40 - 




i 




1000 - 


bOO - 








20 - 






L 


bOO - 




400 ■ 
300 ■ 







3-8 GeV 


1 


Li 


i 



0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 




0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 0.6 0.8 



0.0 0.2 0.4 0.6 0.8 1.0 



0.0 0.2 0.4 



0.6 0.8 1.0 



Figure 3. The energy dependent light curves (histogram) with 
1 - ffl = 0.05, h 1 /h 2 = 0.927 and / = 0.2. The solid lines arc the 
observed light curves from Fermi-LAT (Abdo et al., 2010b). The 
Y-axis is the counts of the observed light curve. 



thicknesses of the main acceleration region to the whole gap 
thickness, /11//12, and the particle number density in the 
main acceleration region, 1 — g\. 



3.2.1 Fractional gap thickness 

Firstly, we consider the azimuthal structure of fractional gap 
thickness and compare the calculated light curves with the 
observations. 

By the definition of the gap fraction of equation 
/ oc — , we may choose the form of the azimuthal distribu- 
tion of / as, 

C 



f(<M 



(22) 



polar cap radius, and the factor 0.18 is chosen by fitting the 
phase averaged-spectrum. 

The pulse profiles using the fractional thickness / de- 
scribed by equation (|22[) are shown in Figure [4] where 
<7i = 0.05 and /ii/ft-2 = 0.927 are the same with those of 
Figure By comparing Figure and 2] we can find the 
effects of azimuthal distribution of / on the pulse profiled. 
For example, the azimuthal structure of the fractional gap 
thickness produces more bridge emissions as well as a third- 
peak-like feature at the phase ~ 0.3 of the pulse profile (in 
particular for higher energy bands) in Figure [4] 

Figure [S] shows the polar cap radius r p (solid- line) and 
the resultant fractional gap thickness / (dashed-line) as a 
function of the polar angle. The polar angle=0° is defined 



Figure 4. The energy dependent light curves (histogram) with 
the distribution of / that provided by the three dimensional outer 
gap. 1 — ji and /11//12 are chosen as 0.05 and 0.927 respectively. 
The solid lines are the observed light curves from Fermi-LAT 
(Abdo ct al., 2010b). 



at the right-hand side of the north pole's magnetic axis in 
the plane including the rotation axis and the magnetic axis, 
where the rotation axis is at the left-hand side of magnetic 
axis in the north pole. Hence the polar angle=180° is defined 
at the left-hand side of the magnetic axis. The angle between 
the magnetic axis and the rotation axis is chosen to be < 90° , 
therefore the outgoing radiation from the south pole can be 
detected with viewing angle < 90°. In the south pole, the 
emissions on the magnetic field lines emerging from polar 
angle larger (or smaller) than 180° produces the first (or 
second) peak. As shown in Figure(5] the fractional thickness 
/ becomes maximum around the polar angle of 200° , where 
r p (4> P ) becomes minimum. Because larger fraction thickness 
/ produces a stronger E\\ as equation (10) indicates, the 
emissions on the magnetic field lines emerging from ~ 200° 
is more stronger than those from ~ 160° and from ~ 240°. 
Consequently, third-peak-like structure is formed in the light 
curves of Figure [4] 



3.2.2 Gap thickness of the two layers 

In the electrodynamic point of view, the thickness of the 
screening region of the gap, hi — hi, will be determined 
by the photo-photon pair-creation rate between the 7-rays 
and X-rays from the stellar surface. It is expected that as 
the null charge surface is closer to the stellar surface, the 
number density of the X-ray photons increases in the gap 
and therefore the screening region becomes thinner. This 
indicates that the ratio of the thicknesses of the main accel- 
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Figure 5. The distributions of the polar cap radius r p (solid 
line) and the resultant / (dashed line) as a function of Polar 
Angles, with the vacuum dipole field. The definitions of the polar 
angle=0° and 180° have been made in the text. 



eration region to the whole gap thickness, /11//12, is closer 
to unity as the null charge surface approaches to the stellar 
surface. Because the radial distance to the null charge sur- 
face is a function of the azimuthal coordinate, it is expected 
that the ratio /11//12 varies in the azimuthal direction. To 
simulate above dependency, let's assume the formula of the 
ratio as 



B. 



l/r„ M »(0p)-l/C° 

1 /„miri 1 / r max 

L l ' null V ' null 



(23) 



where the factor Bi represents the ratio at the polar an- 
gle, at which the radial distance to the null charge surface 
becomes maximum, r™^, and Si + B2 corresponds to the 
ratio at azimuthal angle, at which the radial distance to the 
null charge surface becomes minimum, r™;™. By fitting the 
spectral shape, we obtain B\ = 0.89 and B2 = 0.09. 

To see the effects of the azimuthal dependency of the 
thickness ratio, /11//12, on the pulse profiles, Figure [5] shows 
the simulated pulse profiles with equation (|23p and the fixed 
fractional thickness of / = 0.2 and the number density of 
1 — <?i = 0.05, which are the same with those of Figure [3] 
By comparing Figures and [5] one can see that the distri- 
bution of the /11//12 does not affect much the light curves 
above energy E > lGeV, but it makes a bump in the bridge 
(~ 0.2 — 0.3 pulse phase) of the light curve around lOOMeV. 
We can see those light curves below E < 1 GeV in Fig- 
ure [IJ] are more consistent with the observation than those 
of Figure [3] Wang et al. (2010) show that the ratio of the 
thicknesses hi/h,2 affects the emissivity in the screening re- 
gion. Because the particles in the screening regions, where 
the accelerating electric field is small, mainly produces the 
curvature photons of the energy less than E ^ 1 GeV, the 
azimuthal distributions of the thickness ratio affects the light 
curves below E 1 GeV. 

Figure [7] shows the distributions of the radial distance 
to the null charge surface r nu ii (solid-line) for the inclina- 
tion angle a — 57° and the corresponding /11//12 (dotted- 
line) expressed by equation (|23[1 . It shows that, at the po- 
lar angle~170° that corresponds to pulse phase~0.4, the 
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Figure 6. The energy dependent light curves (histogram) with 
the distributions of /11//12 that provided by the three dimensional 
outer gap. 1—gi and / are chosen as 0.05 and 0.2 respectively. The 
solid lines are the observed light curves from Fermi-LAT (Abdo 
et al., 2010b). 
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Figure 7. The distributions of the distance of the null charge 
surface from stellar surface (solid line) and the resultant /11//12 
(dashed line). 



screening region is thinnest, which means the radiations be- 
low lGeV at the pulse phase ~ 0.4 are weaker than those 
of other pulse phases. Consequently, the shape of the bridge 
emissions of the light curves becomes more consistent with 
the observations and stands out the third-peak-like struc- 
ture. 
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3.2.3 Number density 

The azimuthal distribution of the particles will be more 
complicated, because it is related with (1) the motion of 
the particles under the Lorentz force and (2) the propagat- 
ing direction of the curvature photons, which are eventually 
converted into the pairs in the gap. For example, the non-co- 
rotational electric field perpendicular to the magnetic field 
lines causes E x B drift motion of the particles. Moreover, 
because of the flight time of the photons and of the rotation 
of the magnetosphere, the magnetic field line, on which the 
pairs are created by the pair-creation process, is different 
from the field line, where the photons were emitted by the 
curvature radiation process. 

With the two-dimensional two-layer model, Wang et al. 
(2010) found that the phase-averaged spectra for most of the 
7-ray pulsars can be reproduced by the averaged charge den- 
sity of po = [hipi + (fta — h\)p2\/h2 ~ 0.5. We want to point 
out that because the actual E±((f> p ) is different at different 
4> p , particles in different <^> p -cell may drift into other <f> p -cel\. 
Therefore we assume that the averaged charge density has 
a distribution in the azimuthal direction due to the effect of 
E± x B drift motion. Without the effect of drifting, the num- 
ber of particle in r/> p -cell is No(<j> p ) oc f(4> p )pb(4>p)- With the 
drifting effect, the real number of particle in <f> p -cel\ (N((j>p)) 
may come from (cj> p + A<^> p )-cell, i.e. N(<p p ) = No(<j> p + Acj> p ). 
Since the average charge density p(4> P ) = N (<f) p ) / f (4> p ) , we 
obtain 

+ A0 P ) 

P(<Pp) = Po jT-j—s , (24) 

JVPp) 

where po = 0.5 is the averaged charge density without drift 
motion. It may be expected that displacement of the parti- 
cles due to the drift motion becomes more important on the 
magnetic field line that has a smaller radial distance to the 
null charge surface, because the particles run longer distance 
in the outer gap, which is extending the null charge surface 
to the light cylinder. To take into account such effect, we 
may assume the formula of displacement as 

,l/r„„H(^)-l/C„1f 



A0 P = F- 



(25) 



where F is a fitting parameter, which is chosen as -28°. We 
take A(p p = for r™^, because the null charge surface 
is close to the light cylinder, therefore particles from other 
cells cannot drift in. The negative value (that is, counter- 
direction of the rotation) of the displacement may be ex- 
pected because (1) more particles are concentrated around 
the upper part of the gap and (2) the direction of the non- 
corotational electric field E± is negative value at upper part 
of the gap, indicating the direction of the drift motion is 
counter direction of the rotation. 

Using the relationship between gi, g-2 and /11//12 given 
by equation ([9]) and p = \h\p\ + (fi2 — h\)p2]/h,2, the dis- 
tribution of 1 — gi can be obtained. Figure [5] is the energy 
dependent light curves calculated using the azimuthal distri- 
bution of 1 — gi, where /11//12 =0.927 and f—0.2 are the same 
with those of Figure [3] Compared with Figure [31 it shows 
that the distribution of 1 — gi makes a third peak at the 
bridge region (~ 0.3 pulse phase) when the energy> lGeV. 
This is due to the distribution of the number density in 
the main acceleration region, where the photons of E >1 
GeV are emitted. Figure [9] shows the average charge density 
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Figure 8. The energy dependent light curves (histogram) with 
the distributions of 1 — g\ that provided by the three dimensional 
outer gap. /11//12 and / are chosen as 0.927 and 0.2 respectively. 
The solid lines are the observed light curves from Fermi-LAT 
(Abdo ct al., 2010b). 
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Figure 9. The distributions of the average charge density of the 
two layers (filled-circle) and the resultant 1 — gi (circle) of the 
field fields of different Polar Angles. 



(filled-circle) and the resultant charge density in the main 
acceleration region (circle). The average charge density and 
therefore the number density in the main acceleration re- 
gion acquire the maximum vale at polar angle around 240° , 
which corresponds to the pulse phase around 0.3. 
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Figure 10. The energy dependent light curves (histogram) with 
the distributions of 1 — gi, /11//12 and /, provided by the three di- 
mensional outer gap. The solid lines are the observed light curves 
from Fermi-LAT (Abdo et al., 2010b). 




0.0 0.2 0.4 0.6 0.8 1.0 

Pulse Phase 



Figure 11. Intensity map in the pulse phase and energy plane. 
The darkness represents the percentage of the number of the pho- 
tons of certain interval of pulse phase in the total number of pho- 
tons of certain interval of energy. 



3.3 Combined pulsed profile 

We have discussed the effects of the distributions of/, /11//12 
and 1 — pi on the shape of the energy dependent light curves. 
As discussed in section 13.2.21 the distribution of thickness 
ratio /11//12 mainly affects the light curve of E ^ 1 GeV 
energy bands, and it makes a bump at pulse phase around 
0.25 in the light curve below 1 GeV. The shape of bridge 



emission in the light curves above 1 GeV are mainly affected 
by the azimuthal distributions of the fractional thickness / 
and of the number density 1 — pi. We demonstrated the 
third-peak-like structure can be produced by the azimuthal 
structure of /, /11//12 and 1 — pi. 

Figure [10] shows the pulsed profiles calculated by tak- 
ing into account the azimuthal distributions for all /, h\jhi 
and 1 — pi. We can see in Figure [TOl that there is a bump at 
the trailing part of the first peak in the pulse profile below 
1 GeV, and it shifts towards the second peak in higher en- 
ergy bands. Figure[TT]is the intensity map in the pulse phase 
and the energy plane. The darkness represents the number 
of the photons at the certain interval of the pulse phase. It 
shows that the third peak in the bridge region shifts towards 
the second peak as increases of the photon energy. 

The reason for the shift of the third peak is explained 
as follows. In the light curves in the energy bands lower than 
1 GeV, the distributions of the thickness ratio hi/h,2 makes 
a third-peak structure at ~ 0.2 pulse phase. In the energy 
bands higher than 1 GeV, on the other hand, the azimuthal 
distribution of the fractional thickness, /, and the number 
density, 1 — pi, is more important, and they produce the 
third peak at ~ 0.3 — 0.35 pulse phase. Consequently, the 
differences in the standing phases of the third peak due to 
the distributions of /11//12 and / (or 1 — pi) produce the shift 
of the third peak with the photon energy. 

Here the distribution of the polar radius r p , position of 
null charge surface r nu a, which determine the distribution 
of 1 — pi, /ti//i2 and /, depend on the inclination angle a. 
For instance, if the inclination angle is small, such as 10°, 
the shape of the polar cap is close to a circle, then the r p 
of different polar angles are nearly a same value, the third 
peak made by it is not obvious. On the other hand, if the 
inclination angle is getting larger, the third peak made by 
the distribution of / is getting more obvious. In order to 
provide appropriate distributions of the three parameters, 
we use inclination angle a = 57° and viewing angle /3 = 80° . 

Fermi 7-ray telescope have measured the 7-ray emis- 
sions from about 60 pulsars. However, the Vela pulsar is 
the only one showing the shifting third peak with energy. 
This may be because the Vela pulsar is relatively young and 
its gap thickness is smaller than other mature pulsars. For 
mature pulsars (e.g. the Geminga), the gap thickness will 
be much thicker (/ ~ 0.5) than that (/ ~ 0.2) of the Vela 
pulsar. In such a thick gap, the emissions from the different 
regions in the gap may smear out the effects of azimuthal 
structure of the gap. For the Crab pulsar, the gap thickness 
/ ~ 0.1 is smaller than that of the Vela pulsar. However, 
the contributions of the inverse-Compton scattering from 
the secondary pairs, which are produced outside of the gap, 
dominates the emissions in 7-ray energy bands (Tang et al. 
2008) . This may be the reason that no third peak is observed 
in the light curves of the Crab pulsar. 



3.4 Spectra 

The spectral characteristics are summarized in Figures 1121 
1151 Figures [12] and [13] compare the calculated (solid lines) 
and observed (filled-circles) phase-averaged and phase- 
resolved spectra, respectively. By fitting calculated spectra 
with power law plus exponential cut-off form, Figures [14] 
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Figure 12. The phase averaged spectrum with the distributions 
of 1 — 5i , h\/h,2 and /, comparing with the observed data (circle) 
from Fermi-LAT (Abdo et al., 2010b) 



and [15] show the cut-off energy and the photon index, re- 
spectively, as a function of the pulse phase. 

As we can see from Figure 1141 the present model can 
explain the observed cut-off energy as a function of the pulse 
phase. This can be explained as the effect of the distribution 
of the /, which decides the energies of the photons in the gap 
by determining the value of the accelerating electric field. 
And the cut-off energy of the spectrum is mainly determined 
by the emission from the main acceleration region. From 
Figure [5] in the observable region from around 0° to 260° , 
there are three peaks of /, which correspond to the three 
peaks in Figure 1141 And the pulse phase of the middle peak 
in Figure [14] is consistent with that of the third peak in the 
light curve made by the distribution of /. These two facts 
indicate that the distribution of / makes the calculated cut- 
off energies of the phase-resolved spectra consistent with the 
observed one. 

Because the photon index is much more sensitive to the 
shape of the spectrum than the cut-off energy, the Figure [Pol 
is not so consistent with the observation, where the photon 
indices between pulse phase of 0.3 to 0.5 are different with 
the observed ones. As the spectrum of 0.315 < tp < 0.324 in 
Figure Q2] shows, the present model predicts a larger photon 
index at the pulse phase of 0.3-0.4 than the observations. For 
the pulse phase=0.4 to 0.5, as we can see from the energy 
dependent light curves (Figure HT)]! . the calculated photons 
are more than those of the observation. Consequently, the 
calculated curvature radiation from the main acceleration 
region is too strong, leading the calculated spectrum is hard 
with a photon index ~ 1. For the other pulse phases, as 
shown by the Figure [131 the flux and the shape of the spectra 
are close to the observed ones, and the photon indices of 
these pulse phases can also explain the observed ones. 



4 SUMMARY 

In this paper, a three dimensional outer gap model is built 
to explain the energy dependent light curves of the Vela 
pulsar observed by Fermi LAT (Abdo et al., 2010b). In 





Figure 13. Four phase-resolved spectra provided by the more 
realistic model that using the distributions of 1 — gi, h\/h,2 and 
/, comparing with the observed data (circle) (Abdo et al., 2010b). 
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Figure 14. The cutoff energies of the phase-resolved spectra of 
different pulse phases given by the more realistic model, compar- 
ing with the observed data (circle) (Abdo et al., 2010b). 



the present two-layer model, the acceleration and emission 
region in the gap consist two parts, that is, the main accel- 
eration region and the screening region. In the main accel- 
eration region, the charge density is much lower than the 
Goldreich- Julian charge density and a strong electric field 
accelerates the particles to emit GeV photons via the cur- 
vature radiation process. In the screening region, the charge 
density is higher than the Goldreich- Julian charge density 
to screen out the accelerating electric field. We extend our 
two-dimensional two-layer model in Wang et al. (2010) into 
a three-dimensional one with the azimuthal structure of the 
fractional gap thickness (/), ratio of the thicknesses of the 
primary and whole region (/11//12), and the number den- 
sity in the main acceleration region (1 — gi). Using constant 
1 — gi, hi/Ii2 and /, although the calculated light curves 
can qualitatively explain the observations of the Vela pul- 
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Figure 15. The photon indices of the phase- resolved spectra of 
different pulse phases given by the more realistic model, compar- 
ing with the observed data (circle) (Abdo et al., 2010b). 
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sar, the energy dependent light curves can not have the third 
peak. Therefore, we considered the possible azimuthal dis- 
tributions of these three parameters. We found that the dis- 
tributions of 1 — gi and / make third-peak-like structure in 
the bridge region of light curve of above 1 GeV, while the 
distribution of h\ jhi makes a bump in the bridge region of 
the light curves below 1 GeV. We find that the phases of the 
third peaks caused by the azimuthal distributions of /ii/fe, 
1 — gi and / are different from each other. Consequently, the 
difference in the phases produces the shift of the combined 
third peak with the photon energy. We also showed that the 
present model can reproduce the distribution of the cut-off 
energy for each rotation phase. 
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